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ABSTRACT 



^ ■ A model-independent reconstruction of mechanical profiles (density, pressure) 

of the solar interior is outlined using the adiabatic sound speed and bouyancy 
■^j- . frequency profiles. These can be inferred from helioseismology if both p- and 

o 



g-mode frequencies are measured. A simulated reconstruction is presented using 



Q\ . a solar model bouyancy frequency and available sound speed data. 

ON 



Subject headings: Sun: interior — Sun: oscillations 



O 

?-h ' To appear in the Astrophysical Journal: 10 September 2000 

Standard solar models (SSMs) are now very accurate, exemplified by the work of Bahcall 
et al. (Bahcall et al. 1998; Pinsonneault 1998) and Turck-Chieze et al. (Brun et al. 1998). 
^ ■ But confidence in predicted solar neutrino fluxes is enhanced by checks of the solar core 
depending on only a few general assumptions, independent of detailed models. This can 
be done with helioseismology, although the g-modes are needed. Serious claims of g-modes 
observations have been put forward (Hill & Gu 1990; Thomson et al. 1995) but remain 
controversial (Guenther & Sills 1995; Lou 1996). Nonetheless, the g-modes will undoubtedly 
be observed at some point, despite their attenuation by the convective zone. 

We assume negligible radiation pressure P 7 <C P gas and no significant departures from 
an ideal gas (P = P gas = p^tT//i). 

T = d\nP/d\np , r ad = (d In P/«9p) ad , (1) 
V = rflnT/dlnP , V ad = (<91nT/dlnP) ad . 



1 Contribution to the INT Low- Energy Neutrino Physics Workshop, Univ. Washington, Seattle (June- 
September 1999). 
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r a d = 5/3, V a d = 2/5 for a monatomic gas. The Brunt- Vaisala (bouyancy) frequency N(r) : 



with V M = dlnp/dlnP, is real for propagating g-waves, outside any convective zones. The 
base of the solar convective zone r B cz = (O.7O9)i? (Gough et al. 1996). The adiabatic 
squared sound speed c 2 d = r a( j • Pj p. 

The adiabatic seismology formalism here follows Unno et al. (1989), with some differ- 
ences in notation. The seismic modes are labeled by radial index n (number of radial nodes) 
and angular indices / and m; frequencies degenerate in the last if rotation is ignored. The 
spectrum exhibits "middle" / (fundamental or n — 0, I > 1) modes and the p- (g-) modes 
rising above (falling below) the / mode in v for n > 0, / > 0. For large n, the p- and g- modes 
are concentrated near the surface and the center, respectively. Their frequencies v v and u g 
are 



for large n. (/ cannot be large in v p .) A large set of measured frequencies allow reconstruction 
of the c ad (r) and N(r) profiles. The Bahcall-Pinsonneault 1998 (BP98) SSM is used here 
when specific numerical results are needed. 



Present helioseismic data rule out the T = 5/3 convective polytrope deep into the 
core (Basu et al. 1997; Christensen-Dalsgaard 1997). Hydrostatic equilibrium is tested by 
correct prediction of the v g (n,l), n> 1. Success here would limit any remaining possibility 
of convection at the center (otherwise N 2 < 0), although the Sun is not far from this state 
(M ~ 1.1M in standard calculations; see Kippenhahn & Weigert 1990). An equivalent test 
of hydrostatic equilibrium at the center is provided by the sound speed slope, measurable in 
principle with j9-modes alone: 



N 2 = f • \^ - = (P9 2 /P) ■ [V ad - V + V,] 

*P Li J-ad J 



Ap 1 = -dlnP/dr 



(2) 





1. Tests of Hydrostatic Equilibrium 




1 - 




(4) 



Assume dT^/dr = in the core and a non-singular mass distribution, so that 



g = Gm/r 2 = (4irG/3)p c r — >■ as r — > 



(5) 



Thus dc 2 Ad (ti)/dr must vanish. 
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The measured sound speed slope places a limit on an unperturbed state with non-zero 
velocity field v , taken here as radial for simplicity. This test checks for the presence of 
convective motion that does not efficiently transfer heat. Such convection would not directly 
affect T and the Pj p profile and thus is not ruled out in general by current p-mode helioseismic 
data. Hydrostatic equilibrium is modified to 



dP dv 
dr dr 



(6) 



ignoring spherical-symmetry-breaking rotation and magnetic effects, also calculable. Assume 
again dT ad /dr = in the core and g(0) = 0. Equation (4) no longer vanishes at the center 
and instead yields 



del, 



(0) 



dr 



= -r 



ad 



v (0) 



dv o (0) 
dr 



1 



r(o)J 



Equivalently, in terms of g- modes, v g {n, I) is corrected (to linear order in vq) by the "expecta- 
tion value" in the (n, /) mode of k r v /2n, where k r = radial wavenumber. The related length 
scale ~ i? /few for low-n modes; but for high-n g-modes concentrated near the center, it is 
the central scale R : R% = 3P c /2irGp1, with Rq/R® = 0.12 for the present SSM (Kennedy 
& Bludman 1999). For large n, 



n 



'Bn 3 p c G 



RqQq 



1 



r(o) 



1 

Tad 



, n 9 = / - N(r) 
Jo r 



(7) 



A conservative relative observational error for v g is 0.003, while the relative theoretical 
uncertainty ~ 0.001 (Hill & Gu 1990; Guenther & Sills 1995; Bahcall et al. 1998). Thus a 
direct bound on the core vq < 2%v g {n, l)/k r ~ (7.5 x 10 7 cm sec" - v )^l(l + l)/n\ or 1.2xl0 5 
cm sec -1 for I — 1, n — 30, is feasible, independent of azimuthal splitting by m. 



Mechanical Profiles 



The N(r) and c a d(^) profiles are sufficient to reconstruct uniquely the mechanical (hy- 
drostatic) profiles. N(r) and dcgd(r)/dr are proportional, either one controlling the v g (n, I) : 

dc ad (r) N(r)/2 r. 1 



dr 



Vr ad /r(r) - 1 



1 - 



r(r)J 



ad 



(8) 



Gravity modes are a result of the inhomogeneity of c a( j(r). Assuming dT ad /dr — 0, F satisfies 
a quadratic equation, 

2 



;i + p)r 2 -(2 + p-r ad )r + i = o , v( r ) = 



2 • dc &d (r)/dr 
r ad • N(r) 



(9) 



-4- 



with one physical root, 

r( ) = 2 + V{r) • T ad + y/[2 + V(r) • T ad ]^ - 4[1 + V(r)] 
{ 1 2[l + £>(r)] ' 1 ' 

Thus the T(r) profile is reconstructible with N(r) and a numerical derivative of c ad (r). Note 
that T = 1 when dc^/dr and X? = 0. 

The mass and pressure profiles then follow: 

Gm(r) = - (11) 

{ ) r ad • [i - i/r(r)] ' 



2 



AnGp c = 
P(r) = 



6 • dclM/dr 

T ad • [1 - i/r(o)] 

P(r) -Cad( r ) 

Tad 



where p(r) is obtained from m(r). Note that m > 0, because T > 1 if dc 2 &d /dr < and 
T < 1 if dcl d /dr > 0. (The dc^Jdr = and r — > limits are well-defined in the ratios.) The 
present SSM inner core is an example of the T < 1 case. Dimensionless structure in terms 
of homology variables can also be derived (Kennedy & Bludman 1999). 

Direct methods of inversion without a reference SSM, based on p-modes alone and 
including use of the sound speed slope as a diagnostic of the equation of state, are discussed 
in the older literature (Gough 1984; Christensen-Dalsgaard et al. 1985; Gough & Toomre 
1991; Gough et al. 1995) and, since the advent of SOHO, GONG, and BiSON/LOWL, by 
Gough et al. (1996), Basu et al. (1997), and Christensen-Dalsgaard (1997). Since the z/, 
especially at long wavelength, are affected by the whole Sun, errors in the inferred spatial 
structure are strongly correlated. An inversion for N(r), like the inversion for c ad (r), will 
have uncertainties in both N and the position r. The derived profiles (11) will have errors 
compounded from all these sources: c ad (r), N(r), and r, and thus larger than errors in any 
one. 



3. Numerical Test of the Method 

As the SSM evolves, V falls while V M increases from zero, the core becoming more con- 
centrated. Helioseismology is adiabatic and non-evolutionary and cannot give the thermal 
or chemical structure without further assumptions or a full SSM. Many SSM core features 
can be understood via mechanical-thermal core homology without a detailed model (Blud- 
man & Kennedy 1996). Homology applied to the core only is valid for power-law opacity 



- 5 - 



and luminosity generation and does not require a polytrope. A polytrope is not even ap- 
proximately valid for the present SSM core in any case (Kennedy & Bludman 1999). The 
homology in reality is violated somewhat, as the exponents are not constant and the lu- 
minosity is not exclusively produced by one reaction chain. 2 But an accurate thermal and 
chemical reconstruction is not much less complicated and requires no fewer assumptions than 
a full SSM. On the other hand, the mechanical structure of a SSM can be specified by the 
model-independent reconstruction of Section 2, rather than calculated. 

As p and P are in principle inferable from helioseismology, the numerical precision of 
the reconstruction can be estimated. No g-mode measurements are now available, so the 
SSM bouyancy frequency is used here. Combining c a( j(r), dc^r) / dr , and the SSM N(r) and 
assuming T ad = 5/3, the profile V(r) is obtained (9) and hence T(r). 

Figure 1 shows the measured sound speed below the convection zone, the points being 
fixed by helioseismic inversion (Christensen-Dalsgaard 1997). The BP98 SSM T and N 
profiles, respectively, are shown in Figs. 2(a,b). In Figure 3 is displayed the ratio of the 
reconstructed to the SSM T(r). The helioseismic inversion points are separated by steps of 
A(r/R Q ) ~ 0.01, and we should thus expect agreement with the SSM at level of about one 
percent, consistent with Fig. 3. The inner boundary at r/R Q = 0.05 marks the terminus 
of the present helioseismic inversion, the outer boundary the base of the convective zone at 
tbcz/R& — 0.71. Cutting off the fit at both ends introduces artificial errors, and thus the 
first three and last three points are not shown in Fig. 3. 

The agreement of a true reconstruction with the SSM should be better than Fig. 3, 
but the overall uncertainty would be somewhat larger than a percent, reflecting the c a d(r) 
uncertainty; the additional uncertainties due to a measured, rather than a model, N(r) 
profile; and the errors induced by spatial correlation of the helioseismic inversion mesh 
points. 

4. Summary 

Some solar core tests are proposed, independent of model details. Measured p- and g- 
modes are needed to reconstruct c 2 d (r) and N 2 (r) uniquely and independently of a reference 
SSM, with only general physical assumptions: hydrostatic equilibrium, the ideal gas law, 
and known, constant r ad . 



2 If the specific luminosity production is e — e p x T e and the specific opacity is n = n p q T s , the homology 
yields r e ~ s +3*+3g ^ £qKo . ^ - s - 4 . m e-s+\+q-2 ( correc ti n g an error in Bludman & Kennedy 1996). 
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Absence of core convection implies N 2 > in (2) and the existence of v g {n,l) for high 
n (3). A non-singular, static mass distribution implies vanishing of the sound speed slope (4) 
at the center (5). The measured slope and/or high-n v g can set a limit ~ (9(10 5 ) cm sec" 1 on 
a central velocity field. Reconstruction of the solar interior mechanical structure follows (10- 
11) from c a d(r), dc^ij-) / dr , and N{r). A test of the method is presented, using the available, 
measured c a d(V) and a model N{r) profile. 3 

The c a< j and N profiles to the origin are essential for testing the SSM in the inner 
core, the critical region for the solar neutrino problem, where the higher-energy Be and B 
neutrinos are produced and where solar structure becomes strongly non-polytropic (Bahcall 
1988; Kennedy & Bludman 1999). p-mode data on a finer radial mesh and extending to the 
center are presently lacking, but might become available soon, even if g-mode measurements 
remain elusive. 

Truly independent tests of thermal structure and chemical evolution would require com- 
parison with other Sun-like stars, in particular via asteroseismology (Deubner et al. 1998). 
With the appropriate changes, the method outlined here can be used to analyze such stars 
as well. 

The author thanks John Bahcall (Institute for Advanced Study) and Marc Pinsonneault 
of (Ohio State Univ.) for the BP98 model tables and J0rgen Christensen-Dalsgaard (Aarhus 
Univ.) for helioseismic data. 4 This work is partly based on research done in collaboration 
with Sidney Bludman (Univ. Pennsylvania and DESY) and was supported by the Depart- 
ment of Energy under Grant Nos. DE-FG06-90ER40561 at the Institute for Nuclear Theory 
(Univ. Washington) and DE-FG02-97ER41029 at the Institute for Fundamental Theory 
(Univ. Florida), and by the Eppley Foundation for Research. The author is also grateful to 
the Aspen Center for Physics for its hospitality. 
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FIGURES 

Fig. 1. — The p-mode adiabatic sound speed in core of Sun inferred by helioseismic inversion. 
Uncertainties too small to show. (Data courtesy of J. Christensen-Dalsgaard.) 

Fig. 2. — The BP98 present Sun model profiles from center to base of convective zone at 
r/R & = 0.71: (a) stiffness T, with T — > 5/3 at base of convective zone; (b) Brunt- Vaisala 
bouyancy frequency N, which vanishes at center and base of convective zone. (Model tables 
courtesy of M. H. Pinsonneault.) 

Fig. 3. — Profile of ratio of simulated reconstruction of T to BP98 SSM T, extending from 
r/R Q = 0.075 to 0.68. 
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